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Abstract 

We show that, for pertinent values of the fabrication and control parameters, an 
attractive finite coupling fixed point emerges in the phase diagram of a Y-junction 
of superconducting Josephson chains. The new fixed point arises only when the 
dimensionless flux / piercing the central loop of the network equals vr and, thus, 
does not break time-reversal invariance; for / 7^ tt, only the strongly coupled fixed 
point survives as a stable attractive fixed point. Phase slips (instantons) have a 
crucial role in establishing this transition: we show indeed that, at / = vr, a new set 
of instantons -the W-instantons- comes into play to destabilize the strongly coupled 
fixed point. Finally, we provide a detailed account of the Josephson current-phase 
relationship along the arms of the network, near each one of the allowed fixed points. 
Our results evidence remarkable similarities between the phase diagram accessible to 
a Y-junction of superconducting Josephson chains and the one found in the analysis 
of quantum Brownian motion on frustrated planar lattices. 
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1 Introduction 

Networks of fermionic and bosonic quantum systems are now attracting in- 
creased attention, due to their relevance to the engineering of electronic and 

^ Permanent address: Dipartimento di Fisica, Universita di Perugia, and I.N.F.N., 
Sezione di Perugia, Via A. Pascoli, 06123, Perugia, Italy 



Preprint submitted to Elsevier 



21 January 2009 



spintronic nanodevices. Recently, in Ref.[T], the transport properties of a Y- 
j unction composed of three quantum wires enclosing a magnetic flux were 
studied: modeling the wires as Tomonaga-Luttinger liquids (TLL), the au- 
thors of Ref.p^ were able to show the existence of an attractive fixed point, 
characteristic of the network geometry of the circuit. A repulsive finite cou- 
pling fixed point has been found in Ref. [2], in the analysis of Y-junctions of 
one-dimensional Bose liquids. 

Crossed TLL's are the subject of several recent analytical [3j, as well as numer- 
ical [1] papers: these analyses show that, in crossed TLL's, a junction induces 
behaviors similar to those arising from impurities in condensed matter sys- 
tems. In Ref. [5] it has been pointed out that, in crossed spin-1/2 Heisenberg 
chains, novel critical behaviors emerge since, as a result of the crossing, some 
operators turn from irrelevant to marginal, leading to correlation functions 
exhibiting power-law decays with nonuniversal exponents. 

Impurity models have been largely studied, in connection with the Kondo 
models J6j, with magnetic chains [7], and for describing static impurities in 
TLL's [H]. A renormalization group approach to those systems leads, after 
bosonization [5], to the investigation of the phases accessible to pertinent 
boundary sine-Gordon models [8]. Scattering from an impurity often leads 
the boundary coupling strength to scale to the strongly coupled fixed point 
(SFP), which is rather simple since it describes a fully screened spin in the 
Kondo problem or a severed chain in the Kane-Fisher model [IQl ■ A remarkable 
exception is provided by the fixed point attained in the overscreened Kondo 
problem, where an attractive finite coupling fixed point (FFP) emerges in 
the phase diagram [6]. The FFP is usually characterized by novel nontrivial 
universal indices and by specific symmetries. 

Superconducting Josephson devices provide remarkable realizations of quan- 
tum systems with impurities [TTP^ . For superconducting Josephson chains 
with an impurity in the middle [T31ITT] or for SQUID devices [14|12] the phase 
diagram admits only two fixed points: an unstable weakly coupled fixed point 
(WFP), and a stable one at strong coupling. The boundary field theory ap- 
proach developed in Ref. [TT][T^ not only allows for an accurate determination 
of the phases accessible to a superconducting device, but also for a field- 
theoretical treatment of the phase slips (instantons), describing quantum tun- 
neling between degenerate ground-states; furthermore, it helps to evidence 
remarkable analogies with models of quantum Brownian motion on frustrated 
planar lattices [TSPB] . 

In this paper, we show that, for pertinent values of the fabrication and con- 
trol parameters, a FFP emerges in a Y-shaped Josephson junction network 
(YJJN); then, we probe the behavior of the YJJN near this fixed point by 
computing the Josephson current along the arms of the network. The paper 
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is organized as follows: 

In section 2 we provide a Luttinger liquid description of the YJJN and derive 
the boundary effective Hamiltonian describing the network; 

In section 3 we investigate the fixed points accessible to a YJJN for different 
values of the Luttinger parameter g and of the magnetic field threading the 
central loop of the YJJN; 

Section 4 is devoted to the computation of the current-phase relation of the 
Josephson currents along the three arms of the YJJN, with the purpose of 
determining the current's pattern near the fixed points found in section 3. 
There we evidence the remarkably different effects of phase slips near the SFP 
and the FFP; 

In section 5 we argue that -as it happens with other superconducting devices 
[T7] - a YJJN allows to engineer an effective coherent two-level quantum sys- 
tem, whose states are characterized by two different macroscopic current's 
patterns along its arms; 

Section 6 is devoted to our concluding remarks, while the appendices provide 
the necessary background for the derivation presented in the paper. 



2 Effective Hamiltonian of a YJJN 

The y-shaped Josephson junction network we consider is shown in FigiH It 
is made with three finite Josephson junction (JJ) chains ending on one side 
(inner boundary) with a weak link of nominal strength A and on the other side 
(outer boundary) by three bulk superconductors held at phases ipj (j = 1, 2, 3). 
The three chains are connected by the weak links to a circular J J chain C, 
pierced by a dimensionless magnetic flux /. For simplicity, we assume that 
all the junctions have Josephson energies Ej and X <^ Ej. The Hamiltonian 
describing the central region, Hq, is given by 

_ El j2le^\^'-'t^^''+i] + e-^I^o'-<^r''+|]] , (1) 
2 i=i 

where Ec is the charging energy of each grain, N' is the gate voltage applied 
to the i*'^ junction, and 0q ■* {i = l,2,3;i + 3 = i) is the phase of the supercon- 
ducting order parameter at the i-th grain in C. 
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Fig. 1. y-shaped Josephson junction network: all the junctions are equal to each 
other and have nominal Josephson energy Ej, except for the three ones connecting 
the central region to the endpoints of the chain, that have nominal energy A. 

Following a standard procedure [T3|ll|12j . the Hamiltonian in Eq.([l]) can be 

presented as 

He = -i^ E[^f¥ - ^j:{e'Hs?nst\ + e-^i[St'Y[S?]-} , (2) 
i=i ^ i=i 



with H oc E„ [S^'^Y = - N' - I and [S^'Y = e^''^", where nf,^^ is the 
total charge at grain i (measured in units of e*). 

For H > Ej > 0, the eigenstates of Eq.(l2]) are given by 

• A "fully polarized" ground state: 

|0) = I TTT), with energy eo = -^H; 

• A low-energy triplet of states: 

|1,1) = ^[| iTT) + | TiD + l TTi)], with energy ei,i(/) = -f -E^cos(|); 

|1, 2) = ^[| itt) - e-t| Tit) - e^^l TTi)], with energy e.^) = -f - 
Ejcos(^); 

|1, 3) = ^[| ITT) - TIT) - e-'^ TTi)], with energy ei,3(/) = -f - 
Ejcos(^); 

• A high-energy triplet of states: 

|2, 1) = 73[l iiT) + I iTi) + I Tii)], with energy e^M) = ^ ' Ej cos(|); 
|2,2) = ^[| Tii)-e'^| iTi)-e-t| iiT)], with energy e^M) = f - 
Ejcos(^); 

|2,3) = M Tii)-e-^?| iTi)-e^^| iiT)], with energy e3,2(/) = f - 
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Ejcos(^); 



• A high-energy fully-polarized state |3) = | Hi), with energy €3 = 3H. 

We require that C is connected to the three finite chains via a charge tunneling 
Hamiltonian Ht, given by 



T 



-A ^ cos[ 
1=1 



(3) 



Since X/Ej -C 1, one may resort to a Schrieffer- Wolff transformation [18], to 
derive an Hamiltonian Hb describing the effective boundary interaction at the 
inner boundaries of the three chains. To the second order in A, Hb is given by 



3 3 
Hb ~ B{f) J2 e-''t^(^^^ + A{f) J2 e-'^i\''^i'^'' + h.c. 



(4) 



i=l 



i=l 



where B{f) = g ELi (^^^^(7)), and A{f) = § ELi 



i(fc-l) 



<;o-<:i,fc(/) 



A{f) is, in general, a complex number, equal to —Ewe''^ {Ew > 0). Its phase 
7 is related to the magnetic flux by 



tan7 = — 



E 



(-1)" 

k-2,3 ,,(/) 



eo-ei,i(/) 2 ^fc=2,3 eo-€i.fe(/) 



(5) 



for / = 2k7r and {2k + l)7r {k = 0, ±1, ±2, . . .), 7 = 2A:7r/3 and {2k + l)7r/3, 
respectively. 



The Hamiltonian describing the three flnite chains may be written as 



E. 



L/a 



Ho = ^ E E 



1=1,2,3 i=l 
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L/a~l 
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j=l,2,3 j=l 



-Ejcos(0y^-0w ) + 



d 



N 



d 



—I- 



N 



.(6) 



In Eq. ([6]) 0^*'' is the phase of the superconducting order parameter at grain j of 

the zth chain, — i— 777 is the corresponding charge operator; N is proportional 

84')' 

to the gate voltage Vg applied to each grain, while L and a are the length of 
each chain and the lattice spacing, respectively. E^ accounts for the Coulomb 
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repulsion between charges on nearest neighboring junctions. Following the pro- 
cedure detailed in appendix A, Eq.Q may be written in Tomonaga-Luttinger 
(TL)-form [19j as 



1 /(9$, 



dt 



+ v 



dx 



(7) 



where the fields ^j{x) {j = 1, 2, 3) describe the collective plasmon modes of the 
chains, A = E^-f-^^, v = ^^7^^^^^^^^^, and g = 
with Vf = 271 E J sm{ak f) and kf = aiccos{hEc/ Ej) . 



Vf+AnaA[l-cos{2akf)] ' 



Since at the outer boundary the three chains are connected to three bulk 
superconductors at fixed phases ipj, the fields $j must satisfy the Dirichlet 
boundary conditions 



where j = 1, 2, 3 and nj are integers. On the inner boundary, the three chains 
are connected to C result, one should impose here Neumann 

boundary conditions (i.e., = 0). For our following analysis, it is most 

convenient to introduce linear combinations of the plasmon fields, such as 
X{x) = ^E-=i$,(^), Xi{x) = - $2(x)], and X2{x) = + 

$2(3;) — 2$3(x)]. Since Ew is of order }? jEj^ one has that E-^y/ Ej -C 1 and, 
thus, at X = 0, the fields Xi, X2 also satisfy Neumann boundary conditions. Of 
course, at the outer boundary, X11X2 satisfy Dirichlet boundary conditions. 

In the long wavelength limit, the first term on the r.h.s. of Eq.(jl]) may be well 
approximated as 3i?(/) + const Due to Neumann boundary conditions, 
this term only contributes by an irrelevant constant to Eq.(jl]). Eq.(j4]) may be, 
then, usefully presented in the form 

3 

Hb = -2Ew E ■ cos[a,, ■ x(0) + 7] : , (9) 

i=l 



with ai = (1,0), 52 = (-i Y)' "3 = (-i-^)- The colons (: :) in Eq.Q 

denote normal ordering with respect to the vacuum of the bosonic fields XiyX2- 

1 

The effective coupling Ew is given by E^r = (^jj ® E^. Eq.(^ may be regarded 
as the bosonic version of the boundary Hamiltonian describing the central 
region of a F-junction of three quantum wires, introduced in Ref.[l]. As we 
shall see, setting 7 = (2fc + l)7r/3, allows for the emergence of a new attractive 
fixed point also in the phase diagram of the F-junction of superconducting 
Josephson chains. It should be noticed that this fixed point is attractive, since 
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in a superconducting network, bosons are charged; this should be contrasted 
with the situation arising in F-shaped networks of neutral atomic condensates 
[2j, where the FFP is repulsive. 



3 Phase diagram of a YJJN 



In this section, we use the renormalization group approach to investigate the 
phases accessible to a superconducting YJJN. As evidenced in the analysis of 
other superconducting devices [TT][T2] , there is usually a range of values of the 
Luttinger parameter g for which the phase diagram allows for a crossover from 
an unstable WFP to a stable SEP. For a Josephson chain with an impurity 
P^^fTT] and for SQUID devices (HUE], the crossover is driven by the ratio L/ L^, 
where L is the length of the chain (or the diameter of the superconducting 
loop in a SQUID) and is a pertinently defined healing length [11]. Here, we 
shall show that, when 7 = {2k + l)7r/3, a new relevant boundary interaction, 
emerging in a YJJN at strong coupling, destabilizes the SFP: as a result, since 
the WFP is IR unstable, an IR stable attractive FFP emerges in the phase 
diagram. Remarkably, for these values of 7, the phase diagram of a YJJN is 
similar to the one accessible to a bosonic quantum Brownian particle on planar 
frustrated lattices [15], and to spin-1/2 fermions hopping on F-junctions of 
quantum wires [1]. 



3. 1 The weakly coupled fixed point 



Setting E\Y = defines the WFP, where the fields Xi(a;), X2{.x) obey Dirichlet 
boundary conditions at the outer boundary and Neumann boundary condi- 
tions at the inner boundary. As a result, the mode expansion of Xi is given 
by 



Xi{x,t) = 6 + W- ^cos 

\ 9 n 



TT / 1 



X 



ai(n)_g-if(„+i)t,t 
n + \ 



(10) 



with [ai{n), aj{n')] = Sij [n + |j S^+n'-ifi, = /Ui+27mi2, 6 = pvrnig - nuu], 

(/ii,/i2) = ([(^1 - (P2], -^[{^1 - Vz) - (V^i - V52)/2]), with rii^ = Ui -Uj. 

The perturbative renormalization group equations may be derived from the 
partition function, written as a power series in the boundary interaction 
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strength. From Eq.Q, one gets 



-'0 iV=o ei,...,ejv=±l 
13 ^1-- 



N 



j dri j dT2... J dTN {T^f[ : exp[iej a ■ x{rj)] ^ , (11 

•'=1 



with Zo = n^=o[l - q = exp ]), and /? = (A:bT)-i. In Eq.([II]), 

the lattice step a has to be regarded as the short- distance cutoff, (. . .)o denotes 
thermal averages with respect to Zq, and T,- denotes imaginary time ordered 
products. 

The A^-point functions of the vertex operators : exp[iejdkj ■ xi^j)] : are readily 
computed using Wick's theorem for vertex operators [20]. As Pv/L ^ 1, they 
are given by 



N 



{Tr Yl : exp[iejdk^ ■ x(r,)] :)o = exp 



2 ^ 

9 i<j=i 



(12) 



with 



7^(r,r') = In 



e2L — 



(13) 



As a result, at the WFP, one sees that the scaling dimension of the bound- 
ary interaction in Eq.([n]) is given by hw{g) = 1, and that the dimensionless 

coupling strength G{L) = LE^ scales as G{L) ~ L^^s. 

From the operator product expansion (O.P.E.) between vertex operators 



{: exp [idi ■ xir)] :: exp [idj ■ x(r')] :}^,_ 



V(T — T 



exp [-idk ■ x{t)] : , (14) 



with i ^ j k, one gets the second-order renormalization group equations for 
the complex coupling G{L)e'^"' as 



d[G{L)( 



.11] 



L 

9' 



(15) 



which may be usefully presented as 
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^ = [1 - -m) + 2 cos(37)G2(^) (16) 
at g 

^ = -2sin(37)G^(^) (17) 

( £ = In (■^))- Since Eqs. (11611170 are periodic under 7 — > 7 + resulting 
phase diagram of the YJJN will present the same periodicity. Also, the phase 
diagram strongly depends on whether g < 1, or g > 1. Indeed: 

(1) For g < 1, the linear term in Eq.( fT6l) has a negative coefficient and, thus, 
V7, the system is attracted by a fixed point with G^, = 0. Furthermore, 
Eq.( |T7I) shows that the value of 7 at the attractive fixed point is 7=,, = 
2/c7r/3, if {2k - 1)tt/3 < 7(Lo) < (2A; + l)7r, while it is 7* = {2k + ijvr/S 
if 7(Lo) = (2^ + 1)71/3 0. 

(2) For g > I, Eq. (fT6|) has a positive coefficient; as a result, G{i) grows as 
£ increases. Whether is now finite, or 00, depends on the values of g 
and 7(1^0) • 

In the following subsection, we will derive the perturbative RG equations near 
the SEP. We shall see that, for (7 > | (and for any value of 'y{Lo)) , the system 
is attracted by a fixed point with = 00. For 1 < g < j and for 7(L, 



'0) 



{2k + l)7r/3, the SFP becomes unstable since, for 7(^0) = {2k + l)vr/3, a new 
leading boundary perturbation arises at the SFP. As a consequence a stable 
attractive fixed point emerges in the phase diagram at a finite value of G*. It 
is easy to convince oneself that, for 1 < (7 < | and for 7(^0) 7^ {2k + l)7r/3, 
the stable fixed point is still at G* = 00. 



3.2 The strongly coupled fixed point 



The SFP is reached when the running coupling constant G goes to 00. The 
fields Xj{^)i 3 = 1;2, now obey Dirichlet boundary conditions at x = 0. The 
allowed values of Xi(0),X2(0) are determined by the manifold of the minima 
of the effective boundary potential (Eq.([9])). It is easy to see that: 

(1) for {6k — l)7r/3 < 7 < {6k + l)7r/3, the minima lie on the triangular 
sublattice A, defined by (xi(0), ^2(0)) = {2nni2, ;^[27mi3 + irnu]). 

(2) for {6k + l)7r/3 < 7 < {6k + 3)7r/3, the minima lie on the triangular 
sublattice B, given by (xi(0), ^2(0)) = (27™i2 + ^, -^[2iTni3 + Trnis]). 

(3) for {6k + 3)7r/3 < 7 < {6k + 5)7r/3, the minima lie on the triangular 
sublattice C, given by (xi(0), ^2(0)) = (27^x2 - ^, --^[27rni3 + 7rni2]). 



^ ^{Lq) is the value of the phase 7 at the reference length Lq. It should be noticed 
that, if 7(i^o) = (2A; + l)7r/3, 7 does not scale with L. 
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At 7 = (6fc + l)vr/3, 7 = {6k + 3)n/3, 7 = {6k + 5)n/3, the two sublattices 
A and B, B and C, and C and A become degenerate in energy, respectively. 
From Eq.([n]), one sees that, for 7 ~ {6k + l)n/3, the difference in energy 
between the sets of the minima forming the A and B sublattices is given by 



Ew sin 



7 - 



(2fc+l)7r 
3 



Similar expression hold for the difference in energy 
between the sets of the minima forming the B and C sublattices for 7 = 
{6k + 3)71/3, and for the difference in energy between the sets of the minima 
forming the C and A sublattices for 7 = {6k + 5)it/3. 



The Dirichlet boundary conditions at both boundaries are consistent with the 
mode expansions 



Xj{x, t) =^j + J-{ -—Pj - 2^ sm 

V ^ I ^ n^O 



nnx 



a: 



n —i^nvt 



n 



(18) 



with [a\, a^] = 5*'%5m+n,o- 

For {6k — l)7r/3 < 7 < (6A; + l)7r/3, the eigenvalues of the zero-mode operators 
Pj are proportional to the coordinates of the sites of the sublattice A, and are 
given by 

(Pi,P2)a = \J^g 



2 / ni2 
2^ + 73l"'»-'- 



(19) 



for {6k + l)7r/3 < 7 < {6k + 3)7r/3, they are proportional to the coordinates 
of the sites of the sublattice B, and are given by 



{Pi,P2)b = y2g 



ni2 H — 

27T 3. 



/i2 2 / ni2 

2^ + 711"'' 



(20) 



finally, for {6k + 3)7r/3 < 7 < (6A; + 5)7r/3, they are proportional to the 
coordinates of the sites of the sublattice C, and are given by 



{Pl,P2] 



c 



^9 



■ ^f^i 2- 

nio H 

. 271 3. 



/^2 ^ 2 / ni2 

2^ + 711"'' 



(21) 



The eingenstates associated to the above eigenvalues shall be denoted as 
\ni2,nis)i where i = A,B,C. At the degeneracy points, the merging of two 
sublattices of minima implies a merging of the lattices of eigenvalues of the 
zero-mode operators: for instance, for 7 = 7r/3 the set of the allowed eigen- 
values of {Pi,P2) contains both the values {pi,P2)a and {pi,P2)b, for 7 = tt, 
it contains both the values {pi,P2)b and {pi,P2)c, for 7 = 57r/3, it contains 
both the values {pi,P2)c and {pi,P2)a- 
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At the SFP, one may separately compute the contribution of any one of the 
sublattices A, B and C to the total partition function as 



^ exp I -p 



3eZ 



TTvg 

~17 



ni2 + 



In Eq.(l22D, i = A,B, C, tA = 0, = 1, tc = -1, while = n^=i(l - a;"), 
and q has been defined after Eq. (fTT!) . 

If one denotes by i'i,i'2 the fields dual to xi and X2, one may easily write 
their mode expansion as 

(23) 




with [6'q, Pi\ = i6ij, and j = 1, 2. 



For 7 7^ {2k + l)7r/3, the minima of the boundary potential span only one of 
the sublattices A, B and C. In this case, the leading boundary perturbation 
at the inner boundary is given by a linear combination of the dual vertex 



operators V^, V<^, and V^, defined in terms of the dual fields as 



exp 



±t2j-p,-m 



(J = 1,2,3) 



(24) 



with pi = (0,1), P2 = (^,— |), Ps = they describe instanton 

trajectories connecting two sites in one of the triangular sublattices A, B or 
C ("V-instantons"). The two-point correlation function of the dual boundary 
vertices is given by 



TVVT TVVT 

e L - e L 



(25) 



and, thus, the scaling dimension of Vf'{T),j = 1, 2, 3, is given by hs{g) = ^■ 
As a result, the SFP is stable for g > 3/4 and for 7 7^ {2k + l)7r/3. Thus, for 
3/4 < ^ < 1 and for 7 (2A; + l)7r/3, both the WFP, and the SFP are stable 
and, accordingly, the phase diagram allows for a repulsive FFP. For g > 1 
and for 7 7^ {2k + l)7r/3, the SFP is the only IR stable fixed point: the set of 
the allowed eigenvalues of (Pi,P2) depends upon the value of 7, as discussed 
above. Accordingly, the fixed point partition function is given by Tiilp] in 
Eg. (1221) . for a pertinent choice of i. Remarkably, this shows that the SFP is 
time-reversal invariant, even if the "bare" value of 7 breaks this symmetry. 



Lattice A ■ 
Lattice B • 

Lattice C x 




W-instanton V-'^ 



V-instanton 



Fig. 2. Points on the three triangular sublattices A, B and C: at 7 = (2A; + l)7r/3, the 
energies of two sublattices are degenerate and the minima of the boundary potential 
span a honeycomb lattice, whose sites are connected by W-instanton trajectories, 
shorter than the V-instanton trajectories, connecting sites on the same sublattice. 
The black honeycomb is an elementary cell of the lattice of the minima emerging 
at 7 = 7r/3. 

This is not surprising, though, as the symmetry of the system at the IR stable 
fixed-point is usually higher than the symmetry of the microscopic system. 

For 7 = {2k + l)vr/3, the sets of minima belonging to two sublattices have the 
same energy. As a result, the eigenvalues of the zero- mode operators lie all on a 
honeycomb lattice obtained by merging two triangular sublattices, as sketched 
in Fig 12] for 7 = vr/S, at which point the sublattices A and B merge into a 
honeycomb lattice. The leading perturbation near the Dirichlet fixed point 
contains, now, operators representing "shorter" jumps between neighboring 
minima on the honeycomb lattice ("W-instantons"). 

Following Ref. [15j, one may describe these instantons by introducing an isospin 
operator r, acting on a pertinent two-component spinor [I]- As a result, the 
leading boundary perturbation may now be written as 



^ An |-spinor is associated to a minimum lying on sublattice A and a | spinor to 
a minimum lying on sublattice B. 
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(26) 



i=l 
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with Wj{T 



: exp ^idj ■ iP{t) : and ^ Ej — Ey/. Since the boundary- 
interaction contains the isospin operators r, the relevant O.P.E.'s are obtained 
by combining the multiphcation rules for the isospin operators 



± 



(27) 



with the O.P.E.'s of the bosonic vertex operators 



: e 



ttv{t — t') 



_49 



l±-(r-r')a, 



Terms proportional to , which could be generated to second-order in 
^, are suppressed by the condition Sj=i (^j = 0. As a result, higher-order 
contributions to the /3-function of the running coupling strength C, = only 
appears to order C^^. The RG equation for Q is then given by 

i = [l-/^H^7)]C-2C^ . (29) 



For 7 = 7r/3 the scaling dimension of the boundary interaction, hp{g), gets 
renormalized as 

'-^ = -hAs)e . (30) 



For a small enough value of (, the renormalization of hp{g) may be safely ne- 
glected, since it appears only to the third-order in and one may substitute 
hp{g) in Eq. fl29|) with its bare value ^. Thus, the leading boundary pertur- 
bation at the SEP is irrelevant for (7 > 9/4, while it is relevant for g < 9/4. 
As a result, for 7 = |, there is a range of values of g -namely, 1 < (7 < 9/4- 
where neither the WEP, or the SEP, are stable. The flow diagram then implies 
the existence of a EEP in the phase diagram. In Eigj3l the phase diagram is 
sketched for different values of g: because of the periodicity in 7, only the 
stripe < 7 < 27r/3 is drawn. 

The new attractive EEP emerges as a result of the combined effect of the 
design of the YJJN and of the possibility of tuning the frustration parameter 
7 by setting the dimensionless flux / to tt. Since the circular array C can 
have a very small diameter, self-impedance effects should be negligible. Eor 
F-shaped bosonic networks realized with neutral atomic systems [2], the EEP 
is always repulsive, since those systems are insensitive to external magnetic 
fluxes. 
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G 

Fig. 3. Phase diagram of a YJJN for (from top to bottom): g < 1 (the only IR 
stable fixed point is at G* = 0); 1 < g < 9/4 (the IR stable fixed point is either at 
G,, = oo, or at a finite G=k, according to the initial value of 7); 9/4 < g (the only 
fixed point is at G^, = 00). 

4 The Josephson currents 

In this section, we probe the behavior of a YJJN near each one of its fixed 
points by computing the current-phase relation of the Josephson currents along 
the three arms of a YJJN. 

We find that, for any value of g and for 7 7^ |, the current-phase relation is 
the same as the one of a Josephson junction chain with a weak link analyzed 
in Ref. [TT] while, for 7 = |, one finds new and unexpected behaviors. 
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The Josephson currents in the three arms of a Y-shaped JJN are given by 



h = — 



1 d\nZ 1 (91nZ 

+ 



V2 dni 
1 91nZ 
72 dill 
2 e* 91nZ 
?>(3g (9/i2 



V6 

1 91nZ 

V6 5yU2 



(31) 



where Z is the partition function describing the thermodynamical behavior of 
the YJJN, fl are the phase differences introduced in section 2.2, and e* = 2e 
is the charge of a Cooper pair. In the following, we shall provide the explicit 
form of Eg. (1311) near each one of the three accessible fixed points analyzed in 
section 3. 



4-1 The weakly coupled fixed point 



At the WFP, for g <!■, Hb is an irrelevant perturbation and, thus. 



/3 

Z = Tr|e-^^°T,exp i j dr Hb{t) | 





(32) 



may be safely computed using a mean-field approximation. In Eq. fl32|) . Hq = 
^Y.j=i,2 S^iQ^i("~'^ + l)Q;j('^)) and the boundary interaction Hamiltonian 
has been defined in Eq.(l9]). As a result, one gets 



Z fa Zo exp 



dr {Hb{t))^P 



Zq exp 



2(5Ew 51 cos[ai ■ /i + 7] 



i=l 



,(33) 



where (. . .)o denotes the thermal average with Boltzmann weight e , and 
Zo = Tr[e-'^^«] = l/{n^=o[l - From Eqs. (l3T]l33il . one gets 



h = |sm ai ■ /i + 7 - sm as ■ + 7 | 

h = jsm ^2 ■ + 7 - sm «! ■ + 7 | 

, 2e*Ew 



9 



-{sin[a3 ■ /i + 7] - sin[a2 ■ + 7]} 



(34) 
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Eqs.f l34p explicitly show the dependence of the current's patterns along the 
arms of the YJJN on both the phase differences fl and the parameter 7. 



4-2 The strongly coupled fixed point 



In order to compute the Josephson currents across the three arms of a YJJN 
at the SFP, one has now to account for the contribution coming from the 
zero modes. In order to do so, one should use Eqs. (!3T]) . with the appropriate 
expression for the partition function 7ie[fl] given by Eq. (!22|) . The zero modes 
affect the total energy by the amount 



Tcvg 

~17 



ni2 



(35) 



which is a function of ni2,ni3,/i. At very low temperature and at fixed fl, 
one may approximate the free energy (— ^InZ) with the lowest value of the 
energies £^„^2,ni3 [z^]? given in Eq. (15B]) . For the zero mode eigenvalues belonging 
to sublattice A, for instance, the Josephson currents turn out to be given by 



1. (l^ 

V2 \2n 




+ ni2 



_}_( l^j^ 2^13 + ni2 



V6 

2 //i2 2^13 + ni2 \ 

3 1 27r V3 / ■ 




(36) 



Eqs.( !36|) show the usual [TT] sawtooth dependence on the phase difference 
/i, exhibited by the Josephson current at the SFP. As jl varies within one 
periodicity interval, the integers ni2, ^13 change by ±1. For instance, for — | < 

^ + 5 {5/ii <^ 1), the Josephson 



^<i,fromf^ = f^-5 = -4= 
currents undergo an abrupt jump from 



5to ^ 



h 



e*v 



V2{2n)L \2n 3 



e V 



V2(27r)L V 27r 3 



to 



e*v 



V2{2n)L \27i 3 



e V 



V2{27r)L 



'271 3 



V2, 



e V 



67iL 



corresponding to the shift {ni2,ni3) — > (ni2, riis + 1). 
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a) 




b) 



Fig. 4. a): Behavior of the Josephson current across the three chains near the SFP, 
with the phases /Ui, /X2 chosen as in section 4.2. b): Sketch of the Josephson current 
pattern in the YJJN corresponding to the current patterns across ^2 = fJ'2- 

It should be noticed that, for ^ = —k^ the current in arm 1 switches from 
to a finite value, while the current in arm 2 does the opposite. This sug- 
gests that a YJJN may be useful as a switch commuting between two states 
macroscopically distinguishable by the value of the Josephson current across 
the circuit branches. Finally we mention that, as it usually happens in super- 
conducting networks p^fT^ . the V-instantons near the Dirichlet fixed point 
round off the spikes of the sawtooth function describing the Josephson current 
phase relationship. This effect is discussed in detail in appendix B. 

4.3 Instanton effects for 7 = 7r/3 at the finite coupling fixed point 

For 7 = 7r/3, near the FFP, new more dramatic instanton effects take place in a 
YJJN. Indeed, for 7 = vr/S, two triangular sublattices become degenerate and 
shorter instanton paths are allowed. As evidenced in section 3, the operators 
Wj,Wj representing these paths, become relevant for 1 < < |, and drive 
the system away from the Dirichlet point. Here we evidence the remarkable 
effects of W-instantons on the distribution of the Josephson currents in the 
arms of a YJJN. 

For 7 = f , the minima of the boundary interaction lie on the honeycomb 
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lattice depicted in FiglJl Their position is given by 



2g 



^12 \ 

2 y 



(39) 



with e = 0, 1. Setting /ii ~ /x* + 5 = — | — 6, and — ^ < fi2 < 
with I^I/tt ^ 1, the state | t) = |0,0)yi and the state | |) = \0,0)b are 
quasidegenerate (indeed, they become exactly degenerate for 5 = 0). For 
this choice of the quasidegenerate states, the W-instantons are described by 



1^(0) = -2C : cos 



write their contribution to 



. Substituting Eq. flB.6p into Eq. flB.7p allows to 
the partition function as 



dx \ ix + Co + sg{a)a5 — C^r[l — 2hp{g)][ix + Cq + sg{a)a5] 

271 I 



2V{9)- 



[ix + e,f - + C'] - C'r[l - 2hp{g)] Y^[ix + eo + ^g{^)a5f^-^^^ 

7 

- C^r^il - 2h,{g)][{ix + e,f - a^T'^^^-'] , 

where a = g/{QTi), cq = eo(5,/i2) = ^(| + i^ + 0),a; = f^u;, and sg((T) = 1 
if cr =1, = -1, if a =[. 

To compute Eq. fl40l) is quite a formidable task: however, an approximate com- 
putation can be carried out, for = | — e (e <^ 1), near the FFP ^=1. = (1 — ^)2. 
Indeed, since = ~ ^ 1, neglecting 0(C*)-terms in Eq. (HO|) . leads to 



Z[5,/i2] ~ exp 



cosh 



cosh 



(41) 



from which, for jSv/L ^ 1, one gets 



e f 



71 



e*v 2 



2^,2 



^252 + 3^2 
^252 + 3^2 



/^2 



(42) 



Eqs.f l42p yield the current-phase relations near the attractive FFP shown in 
Figini The typical sawtooth behavior of the Josephson current-phase relation 
is now associated to a stable attractive FFP in the phase diagram. 
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a) 




Fig. 5. a): Behavior of the Josephson current across the chains 1 and 2 near the FFP, 
with the phases /ii,/X2 chosen as discussed in section 4.3. b): Sketch of the Joseph- 
son current pattern in the YJJN associated to the macroscopic states |c(=T)i)) 
introduced in section 4.3. 



Usually, in superconducting systems, such as SQUIDs [1^1[T2] and Josephson 
chains with localized impurities p^lITT] , either the smoothening of the spikes of 
the sawtooth function describing the Josephson current-phase relationship at 
strong coupling is a perturbative effect, or the SFP is unstable, since quantum 
fluctuations drive the system to the WFP. At variance, for a YJJN at the 
FFP, the smoothening of the spikes of the Josephson current due (now) to the 
W-instantons is a nonperturbative effect and the FFP is a stable attractive 
fixed point in the phase diagram. Since a sawtooth behavior of the Josephson 
current is usually associated [T^lllll2j to the emergence of a macroscopically 
quantum coherent two-level system in the superconducting device |T7], one 
may safely expect that an effective macroscopic two-level quantum system 
-this time robust against quantum fluctuations- may emerge in a YJJN, as 
well. This issue has been addressed in Ref.pl] and will be revisited in the 
next section. 
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5 A quantum two-level system emerging in a YJJN 

In this section, we revisit the arguments given in Ref.[21] to show that an ef- 
fective quantum two-level system with frustrated decoherence [22] may emerge 
from a YJJN near the FFP. 

As evidenced in appendix A, the energy of the long-wavelength excitations of 
the YJJN is given by 

E=^[P\' + E' , (43) 

where p = (^1,^2) is the eigenvalue of the zero-mode operators P = {Pi,P2), 
introduced in section 3, and E' accounts for the energy of the plasmon modes 
described by the TLL Hamiltonian given by Eq.(l7]). Plugging Eqs. (ll9ll20p into 
Eq.( H3i) yields the explicit dependence of the energy on the minima on the 
phases ipj of the three bulk superconductors. 

It is easy to see that, for any value of 7 and for all possible values of the 
Luttinger parameter g, it is always possible, for a finite YJJN, to choose the 
phase differences /ii and fi2 to obtain two low-energy quasidegenerate states, 
well separated from the rest of the spectrum. Slightly generalizing the notation 
introduced in Section 4, we still denote the two quasidegenerate states by | t) 
and I I) and observe that, for 7 7^ {2k + l)vr/3, both states belong to the same 
triangular sublattice (A, B or C), while, for 7 = {2k + l)n/3, they belong - as 
in Section 4- to two different sublattices. 

The dynamics of the two states | t) and | |) interacting with the plasmon 
modes residing on the three chains of the Y-j unction may be written as 

H2 = [E^{lI)-E^{fI)]a^-YO{0)a- -YO^{0)a^ , (44) 

where O is one of the vertex operators Vj (if the states lie on the same 
triangular sublattice), or Wj (if the states lie on the honeycomb lattice ob- 
tained by merging two triangular sublattices), and = Z]cr=ta |sg(cr)|o-)(cr|, 
(T"*" = I t)(i I) o"~ = I i)(T I- In Eq. fH4|) . E^{fl) and Ei{jl) are the energies as- 
sociated to the states | ]) and | |), respectively. The cj^-term contributes only 
if the I t) and the | [) states are quasidegenerate, which may be achieved by 
a slight detuning of the phase differences /^i and /i2 by an amount 5 (^ 2ti). 
The terms proportional to Y describe an effective field in the x-direction and, 
at the same time, the coupling between the transverse components of the spin 
and the bath provided by the plasmon modes of the three chains: on one hand 
they determine a y-dependent renormalization of the energies of the effective 
two- level system -the tunnel splitting of the energies of the states |ct)-, on the 
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other hand, they may lead to the formation of an entangled state between 
the two-level system and the bath formed by the plasmon modes in the net- 
work. This latter effect is a main source of decoherence in a two-level system 
interacting with one (or more) baths j22pT] . 

Depending on whether 7 7^ (2/c-|- l)7r/3, or 7 = (2fc-|-l)7r/3 and on the value of 
the Luttinger parameter g the interaction of the system with the bath provided 
by the plasmon modes of the network leads to different coherent behaviors of 
the YJJN [21j. In the following we shall compute the spectral density of the 
two level system near the SFP and the FFP; as pointed out in Ref . [22] . the 
spectral density provides a measure of the amount of entanglement between a 
two level system and the pertinent environmental modes. 



5.1 Spectral density of the two-level system near the strongly- coupled fixed 
point 



As evidenced in section 3, for 1 < g < 9/4 and 7 7^ {2k + l)vr/3, or for 
(7 > 9/4 and V7, the YJJN exhibits an IR stable SFP in its phase diagram. If 
— 7r/3 < 7 < 7r/3, Eq. ([T^ implies that the quasidegenerate states | ]) and | [) 
lie on the triangular sublattice A and that Eq. (H^ may be explicitly written 

as 

H2,sYV = [E^{fi)-E^{fi)]a'-YV{<d)a~ -YV\<d)a+ , (45) 



where V^(0) is the V-instanton vertex operator. 

To compute the spectral density of the two- level system in Eq. psj) one needs to 
evaluate x^(f2)/n vs. fi, where xl(^) is the imaginary part of the transverse 
dynamical spin susceptibility. Since V-instantons are an irrelevant perturba- 
tion, by neglecting higher-order corrections in Y (see appendix B), one gets 

oc 5(fi - 2A(/i)) + 5(fi + 2A(/3)) , (46) 



with A(/I) = -y/[i?|(/I)]2 + F2_ YioTd Eq. (H^ one sees that the spectrum of 
Eq. fHSl) is given by two classical states, with VL = ±A(/2). As pointed out in 
ref. [21] this behavior signals that there is no entanglement between the two 
level quantum system and the plasmon modes. Since Y is irrelevant (i.e., its 
fixed point value is = 0), there is not even tunnel splitting between the two 
degenerate states, no quantum coherence may emerge in this regime [21] 
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FFP 
DFP 
NFP 



Fig. 6. Behavior of the spectral density x±{^)/^ vs. Q, near the WFP, the SFP, and 
the FFP. 

5. 2 Spectral density of the two-level system near the finite- coupling fixed point 

For 1 < g < 9/4 and 7 = 7r/3, the two states | t) and | |) lie on nearest 
neighboring sites on the honeycomb lattice obtained by merging the sublat- 
tices A and B. As evidenced in section 3, short W-instantons are a relevant 
perturbation at the SFP and render the FFP IR stable. Near the FFP, the 
two-level system is described by 

H2,F^P = [E^{lI)-E^{fI)]a'-^W{0)a~-^W\0)a+ , (47) 

where, now, W{0) is a W-instanton operator. The computation of the spec- 
tral density xl(^) is detailed in appendix B using a self-consistent RPA ap- 
proximation. As a result, the spectral density has now two peaks centered at 
±A,(/i) = ±y/[E|(/I)]2 + (C,/L)2, with a finite width oc f (C)^+^^ where C 
is the finite fixed point value of the running coupling constant, determined 
in section 3. The spectral density is plotted in FigEl where we report, for 
completeness, also the spectral density arising near the WFP |21] . 



6 Concluding remarks 

We showed that, for 1 < g < 9/4 and for / = vr, an attractive FFP emerges 
in the phase diagram accessible to a YJJN. The new fixed point does not 
break time-reversal invariance, and it is a stable, attractive fixed point only 
when the dimensionless flux threading the cental loop equals vr; for / 7^ vr, we 
have shown that only the SFP survives as a stable attractive fixed point. Our 
results show remarkable similarities between the phase diagram accessible to 



22 



a YJJN and the one arising in the analysis of the quantum Brownian motion 
on a frustrated triangular lattice p3][T6] . 

Crucial to our analysis is the fact that, at / = vr - for 1 < g < 9/4 - the 
W-instantons become a relevant perturbation and, thus, destabilize the SFP. 
These instantons emerge ultimately as a result of the Y-shaped geometry of 
the network, since they arise when the minima of the boundary potential span 
the honeycomb lattice depicted in Fig|2j Intuitively, they may be regarded 
as the result of the "deconfinement" of the V-instantons in its elementary 
constituents happening- when 1 < g < 9/A- only at / = vr. A Coulomb gas 
approach could be a very helpful tool to further clarify the nature of the FFP. 

We computed the current-phase relations along the arms of the YJJN near 
each one of the allowed fixed points. We evidenced the parameter regions 
where a YJJN may be operated as a Josephson switch and we showed the 
different effects of the instantons on the current pattern near the SFP and the 
FFP. In particular, in a YJJN at the FFP, the smoothening of the spikes of 
the sawtooth dependence of Josephson current on the phase differences /I is a 
nonperturbative effect, due to the attractive nature of this fixed point. 

Finally, we provided additional arguments confirming that, near the FFP, a 
YJJN supports a quantum coherent two-level system with frustrated decoher- 
ence [2T] . 

In order to set a YJJN to be a quantum device either acting as a Josephson 
current switch or modeling an effective two level quantum system one needs, 
first of all, to promote the phase differences /ii,/i2 to control parameters. 
This may be achieved by resorting, for instance, to multipolar magnetic coils 
[23] inserted in external loops connecting the bulk superconductors at the 
outer boundary of the YJJN: indeed, for sufficiently long chains, the localized 
magnetic fields generated by the multipolar magnetic coil may be engineered to 
avoid variations in the flux threading the circular Josephson junction array C. 
Furthermore, when the YJJN has a finite size L, it is easy to convince oneself 
that the FFP is stable against small fluctuations of the flux /, provided that 
v/L is sufficiently big: for instance, if the point 7 = vr/S is displaced by a small 
amount u, v/L needs to be larger than the energy splitting sin(z/) between 
the minima of two triangular sublattices. At variance, when v/L < Ew sin(z/), 
there is a flow towards the SFP and, depending on sgn(z/), the minima of the 
boundary potential lie on either one of the triangular A and B sublattices [21j. 
Finally, today 's technology allows to fabricate superconducting devices with 
values of g ranging from (7 < 1 , to (7 ~ 2 [21] . 

Josephson networks where n finite chains are connected to a central circular 
array C may be analyzed with tools similar to those used in this paper. Of in- 
terest is also the JJ network with n = 4 since it corresponds to the tetrahedral 
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qubit proposed in Ref . [25] . 
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A Tomonaga-Luttinger description of superconducting Josephson 
junction arrays 



Here, we briefly review the derivation of the effective TLL Hamiltonian, de- 
scribing one-dimensional arrays of Josephson junctions. For this purpose, in 
Eq-dS]) one should assume that Ej/ Ec <ti I and N = 2n+l + 2h, with n integer 
and -C 1 [T31ITT] : then, if one defines effective lattice spin- 1/2 operators as 

Pg , Sf = PGe^"^^PG , (A.l) 



G 



d 



N 



with Pg the operator projecting onto the subspace of the charge eigenstates 
with the charge at any site either equal to or to n + 1, one may present 
Eq.© as 



-PG-^chain-PG = 



E 



L/a-l 



spin 



^ i=i 



+ 



E., 



3 {Ejf 
16 E, 



L/a—l L/a 

EQZ QZ TT \ ^ QZ 



(A.2) 



Eq. (lA.2p is the Hamiltonian for an XXZ-chain in an external magnetic field 
H = hEc P^IITT] : to map it onto an effective TLL Hamiltonian, one needs to 
write the spin operators in terms of lattice Jordan- Wigner fermions aj. Upon 
defining the lattice Fourier modes as 




(A.3) 
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-Pc-^chain-PG in Eq. (jX2]) is given by 



Hjw = ^[-Ejcos{ka) - H]alak + 

k 



3 {E. 



J) 



21 L/a-l 



16 



E 

J i=l 



2) .(A.4) 



From Eg. flA.4p . one see that two "band- insulating" phases open up when > 
Ej [T5][TT] . In spin coordinates, they correspond to fully polarized spin phases 
which are the Coulomb blockade insulating phases setting in the chain when 
the gate voltage is tuned far from charge degeneracy point. 



For \H\ < Ej, Eq. flA.4p describes a one-dimensional conductor. By keeping 
only long-wavelength modes around the Fermi points = ±^arccos{H / Ej), 
and by bosonizing Eq. flA.4p one gets the Sine-Gordon Hamiltonian 



— dx 

4:71 J 





1 fd^y fd^' 
V \dt I \dx 



Gu J dx cos[2V2g^{x) + Akfx] , (A.5) 



3 iEj)2 
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with the Luttinger parameter defined in section 2 and Gu (x E^ — 
When g > 1/2, the last term in Eq.l lA.Sp may be neglected, in the thermo- 
dynamic limit and H'' reduces to the Hamiltonian of a spinless TLL. g may 



either be < 1, or > 1, depending on whether A( 
sive TLL), or A < (attractive TLL) [n]. 



E' - 



3 (Ejy 



16 Ec 



> (repul- 



The normal modes of a spinless TLL may be constructed by introducing the 
dual field ip, related to $ by = ^ and 7;^ = f^, and by introducing two 
chiral bosonic fields, (f)R,(f)L' 



V at 
as 



In terms of 0r, 0l, is given by 



= — I dx 

4:71 




+ 



dx 



(A.7) 



The normal mode expansion of (pR^x — vt), (j)L{x + vt) may be written in terms 
of the Fubini-Veneziano chiral fields ESI as 



L ~zl n 
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(PL{x + vt) = qR + ^PL{x + vt) + zY.^^e^'"'^^+^''> , (A.8) 

with 

[qr, Pr] = [ql, PL]=i ; [aR{n), aR^m)] = -[aiin), aiim)] = n6n+m,o , (A.9) 

with all the other commutators vanishing. As a result: 

H' = -riiPRf + (Pl?] + 1- E[««(-^)«i?H + o^Ma^i-n)] .(A.IO) 



To construct the Fock space, one needs to define a vacuum \{j)r,Pl),^) for 
any allowed pair of eigenvalues of the zero-mode operators Pr, P^, and then 
act with creation operators ^/{(n), n) {n < 0) on the states \{pr,Pl),0), 
which obey the conditions 

Pr\{Pr,Pl),0) =pr\{pr,pl),0) , Pl\{pr,Pl),0) = pl\{pr,Pl),0) , 

aRin)\{pR,pL),0) = aLi-n)\{pR,pL),0) = (n > 0) , (A.ll) 



In a system with boundaries, the boundaries conditions may be accounted for 
by means of pertinent relations between the r and the l modes. For instance, 
Neumann boundary conditions at x = 0, that is, = 0, imply 

Pr-Pl = , aR{n) + «L(-n) = , Vn , (A. 12) 

while Dirichlet boundary conditions at x = 0, that is, $(0) = 0, imply 

Pr + Pl = , ttflH - aL(-n) = , Vn . (A. 13) 



B The partition function and the spectral density of the effective 
two-state system 



Here we set up the general formalism needed to include the instanton contri- 
butions to the partition function of the effective two-level system described in 
section 5. In doing so, it is most convenient to write the spin-1/2 operators 
introduced in Eg. (14^ by means of two pairs of fermionic operators, aa,a]^, 
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such that cr^ = | X]cr=tJ, sg(c'")aJ.ao-, 0"+ = aja^. When doing so, the imaginary 
time action of the effective two-level system reads as 



d_ 



+ [a\a^O{T)+h.c]^ + Sf ,(B.l) 



where S^^^ is the Euclidean action for the plasmon field, given by 

/3 L 



E - 

i=l,2,3^" 



dr J dx 





1 /(9$. 



V \ dr 



+ v 



dx 



(B.2) 



while the chemical potential is iujo = [22] ■ 

At low temperature ^ 1), one may approximate the partition function 
of the effective two-level system as 

ZEff^EZ.<x(/?) , (B.3) 



with Z,,,(r) = (a,(r)a;,(0)). 

The diagrams used to compute Zg-o-' (uj) = /q°° dre^^'^'^Z^j^' (r) are schematically 
depicted in Fig JB.lj there, the solid thin line corresponds to the propagator 
of a fermion a^, given by g^\uj) = ^/[iiuj — coq) + E(j(jl)], and the dashed 
line corresponds to the propagator for the F-vertex. The pertinent Dyson's 
equations yielding Zo-o-'(i^) are given by 



(B.4) 



where a is | (|) if cr is | (t), 7a(co') corresponds to the "bubble" diagram in 
Fig jB.ll given by 
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Fig. B.l. Diagrams contributing to the Dyson equations for Z|| and for Zj-f. 

with h{g) = hs{g) = for V-instantons, while h{g) = hp{g) = \g for W- 
instantons. When Luo/v ^ 1, Eq. (]B.5|) may be approximated as 



/o r \ 2/i(g) 

7.H^ — T[l-2h{g)][i{uj-uJo)+EM?'^^'^-' ■ (B.6) 



From Eqs. (lB.4p . one obtains 



z..H = {[#]-iH-yVM}/ 



+ FVH7<.M .(B.7) 



Eq. (]B.7p has been used in section 4 to compute the Josephson currents near 
the FFP. If 1 < 5( < 9/4, for 7 = {2k + l)7r/3, O is the relevant W-instanton 
operator while, for 7 7^ (2A;+l)7r/3, O is the irrelevant V-instanton operator. In 
the latter situation, the (9(y^)-approximation to Eq. (1B.7l) allows to compute 
the smoothening induced by the V-instantons on the sawtooth behavior of 
the Josephson current. Setting, for instance, — | <7<f)2^~i^^~ ""73' 
and — I < ^ < |) the states | ]) and | J.) belong both to sublattice A, and 
O = —2Y : cos 2^1^2(0) :. In this case, Eq. flB.7p may be approximated as 
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The partition function is then given by 



Z[/ii,/i2] ^ 2exp <^ -(3 



girv 



27r/ 



+ 



V 27r 



1 

+ 3 



2tt 



2Txvg 



,(B.9) 



with y = LYr^ L^-hs{g)_ prom Eq. flRQjl . one may derive the Josephson current 
distribution in the three arms of the YJJN 



'/i2 — ^2 



L \2tiJ V 2 



in 



1 

1 + 2 



e 



x/2 V27r/ V6 V 27r 



~ /^2 

27r 



V 27r 



+ 



£^ /2 //X2JV^ 

' L V 3 V 27r 



A^2 ~ /^2 

27r 



+ 



(B.IO) 



The formahsm developed in this appendix allows also to compute the (trans- 
verse part of the) dynamical spin susceptibility of the emerging two-level sys- 
tem, x±{^)- As discussed in section 5, the imaginary part of x±{^)/^ is the 
place to look at, in order to analyze the entanglement between the system and 
the environmental modes. 

The starting point to derive the xx and the yy-components of is the 

computation of XiTi^) ^^"^ Xi^i^)^ thai is, of the Fourier transforms of the 
imaginary time dynamical susceptibilities XiTi^)^ Xi^i'^): respectively 
given by 



Xr(r) = (T.[a+(r)a-(0)]) , xl^(r) = (T.[^-(r)a+(0)]) 



(B.ll) 



The approximate computation of x±^{^) is graphically shown in Fig JB.2b ). 
To lowest order in Y, x±~{^) is computed as a loop defined by the | t)-state 
propagating forward in (imaginary) time, and by the | |)-state propagating 
backward, while xT^i^) is computed in the same way, by just exchanging | 
and |. Accordingly, x±~{^) and x±^{^), are given by 



1(0), 



J ^Z,,(^)Zu(^ + fi), [xl+](°)(^fi) ^ / ^Z,,(^ + fi)Zu(.;),(B.12) 
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where the functions Zo-o-(c<j) have been defined in Eq. (lB.8p . As a result, one 
obtains 



with A(/i) = + F^. A similar formula holds for [x± pro- 

vided one exchanges t with |, and vice versa, in Eq. flB.13p . For 1 < g < 9/4 
and 7 7^ (2fc + 1)it/3, and for (7 > 9/4, one may safely neglect higher order 
corrections in Y to so that Eq. flB.13p provides a reliable estimate of the 

transverse dynamical spin susceptibility. Both the xx and the yy components of 
X±{^) are obtained from Eq. (]B.13p . and from the analogous one for 
Their imaginary part is computed via the replacement Q — > —iQ + 0"'': in 
both cases it is equal to xK^)) given by 

^^^^ (x6{Q-2A{fI)) + 6{Q + 2A0)) . (B.14) 



Eq. flRTil) is the estimate of 

near the SEP, quoted in section 5.1. 

Eor g < I, the instantons provide a relevant perturbation: thus, higher-order 
contribution in Y to Eq. (1B.12l) cannot be neglected. By taking the large-F 



limit of the fully dressed expression for Z„^{uj) derived in Eq. flB.7p . one gets, 
for the imaginary part of the transverse dynamical spin susceptibility 

X±i^) no IT r7r\ , ni3-fg _ \ow.JTr\ _ 0|3-f ffl 



(x[\2E^{il) + nY--3 -\2E^{pi)-nY--^]in . (B.15) 



Eq.f lB.15l) shows that the largest part of the spectral weight is now in the 



region around f2 = 0: this signals the onset of a fully entangled state between 
the two state system and the bath formed by the plasmon modes [2T|22] . 

For 1 < (yf < 9/4 and 7 = {2k + l)7r/3, the behavior of the system is ruled 
by the IR stable FFP. An estimate of may now be done, for instance, 

when (yf = I — e, with e <^ 1: since the FFP is at C* ~ ea, using the effec- 
tive Hamiltonian in Eq. fl47|) . one may resort to the RPA computation of the 
dynamical spin susceptibility, graphically drawn in Fig lB.2b ). to get 

1 



l]-A,(/I)-er[-l-|e](-(]) 
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b) 





+ 




-L RPA _L 



[xW 



-L RPA 



Fig. B.2. Diagrams contributing to the RPA computation of the dynamical spin 
susceptibihty: the dashed hue represents the propagation of the W-instanton. 



+ 



1 + le 



(B.16) 



with A,(/2) = ^[E|(/3)]2 + {C/iy. Computing From Eg. lfRTHl) . one 

sees that, on one hand, the energies of the two-level quantum system are 
renormalized by ( to ±A^,{fI), on the other hand, that the two peaks at the 
renormalized energies now display a finite width oc ^(C*)^^^^? which is the 
result quoted in section 5.2. 
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